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FOSEtfQiS 


Analytical  studies  on  the  optimum  propellant  combinations  for  large 
multi-stage  rockets  have  clearly  indicated  the  superiority  of  liquid  hydrogen 
as  fuel  for  the  upper  stage  propulsion  units. 

Lack  of  experience  and  technical  know* how,  largely  the  absence  of 
efficient  and  sufficiently  reliable  ancillary  equipment,  3Uch  as  tanks,  pumps, 
valves,  control  equipment,  etc.,  inhibited  the  use  of  hydrogen  in  the  first 
generation  of  ELDO  vehicles. 

To  fill  this  gap,  various  research  programmes  were  launched  in  this 
country  about  4  years  ago  at  EHDE,  EPE  and  RAE,  with  the  objective  to  gather 
fundamental  data  and  technological  information  on  various  aspects  connected 
with  the  hydro  gen/oxygen  propulsion  system. 

This  Establishment  undertook  to  carry  out  extensive  heat  transfer  studies 
with  liquid  hydrogen  in  its  liquid  and  near-critical  fluid  phases.  In  order  to 
make  the  results  directly  applicable  to  the  requirements  of  the  design 
engineers,  test  sections  were  developed  having  the  shape  and  sise  of  cooling 
channels  likely  to  be  used  in  a  second-stage  propulsion  unit  of  a  three-stage 
rocket  vehicle. 

Extensive  preparatory  studies  and  surveys  leading  to  the  design  and 
construction  of  a  complex  test  rig  were  undertaken  during  1966  and  1967* 

However,  shortly  before  the  actual  commencement  of  experiments  with  liquid 
hydrogen,  the  change  in  this  country’s  space  research  policy  resulted  in  the 
cancellation  of  almost  all  work  involving  liquid  hydrogen  for  rocket  propulsion, 
and  consequently  in  the  abandonment  of  this  projected  research. 

The  work  involved  the  full-time  effort  of  several  members  of  the  staff  at 
ERDB  for  over  2  years,  and  considerable  material  expense.  It  was  thought 
desirable,  therefore,  to  record  the  status  attained  in  detail,  should  it  he 
decided  to  take  up  this  work  again  at  some  future  date. 

In  the  first  part,  an  exte^ive  survey  of  all  known  published  information 
on  the  subject  of  heat  transfer  to  fluid  hydrogen  at  supercritical  pressures  is 
presented;  in  the  second  part,  the  aim:*  of  the  projected  research  are 
discussed  and  a  comprehensive  description  cf  the  design  details  of  the  apparatus 
and  its  associated  equipment  are  given,  including  the  proposed  operating 
procedures  and  safety  features  of  the  test  rig. 


H.  Ziebland 
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SUMMARY 


The  published  experimental  data  covering  forced  convection 
heal  transfer  to  cryogenic  hydrogen  are  reviewed,  with  special 
attention  to  the  near-critical  regions  of  temperature  and  pressure. 
Data  for  straight  and  curved  tubes,  of  both  circular  and  non-circular 
cross-sections,  are  covered;  also  the  case  of  asymmetric  peripheral 
heat  flux  through  the  walls. 

A  number  of  theoretical  and  semi-empirical  treatment?  of  the 
near-critical,  variable  fluid  property  condition  are  discussed,  and 
their  effectiveness  in  correlation  of  near-critical  heat  transfer  to 
cryogenic  hydrogen  considered. 


Reference:  WAC/1 90/05 


*.  INTRODUCTION 

The  feasibility  of  using  hydrogen  under  supercritical  conditions  as  a 
regenerative  rocket  coolant  has  been  established,  both  by  experiment  and  by 
computation  (l)  based  on  a  knowledge  of  the  physical  properties  over  the  range 
of  likely  operating  conditions,  including  cryogenic  temperatures. 

Experiences  of  hydrogen  rocket  design  in  the  USA  suggest  that  the 
behaviour  of  hydrogen  as  a  coolant  is  more  difficult  to  predict  using  the 
conventional  forced-convection  correlations,  and  allowing  for  the  variations 
of  fluid  properties  rath  temperature  and  pressure,  than  that  of  most  fluid 
coolants. 

Hydrogen  is  an  unusual  regenerative  coolant  in  that  it  is  possible  for 
the  fluid,  while  at  supercritical  pressures,  to  be  heated  from  subcritical 
bulk  temperatures  (since  storage  may  be  at  temperatures  close  to  the  normal 
boiling  point  at  20°K),  through  its  critical  temperature  (33°K)  during  passage 
through  the  coolant  channels. 

At  some  locations  in  the  coolant  passages,  therefore,  the  bulk  fluid 
temperature  will  attain  values  in  the  region  above  the  critical  point,  in 
which  the  physical  properties  of  the  fluid  are  very  strongly  temperature 
dependent  and  the  specific  heat  (0  )  and  Prandtl  Number  show  marked  peaks  on 
individual  isobars  for  pressures  immediately  above  the  critical  pressure. 

It  is  therefore  necessary  to  determine  how  the  convective  heat  transfer 
processes,  and  in  practice  heat  transfer  coefficients,  differ  from  those 
adequately  represented  by  the  well-established  constant  property  correlations; 
and  desirable  to  provide  a  means  of  reliably  estimating  heat  transfer  in  actual 
cooling  channels  of  rocket  engines  designed  to  operate  under  these  conditions. 

As  a  consequence  of  operating  in  the  critical  region  of  temperature  and 
pressure,  large  changes  in  fluid  density  will  occur  along  some  sections  of  the 
coolant  tribes,  and  these  must  also  be  estimated  reliably  for  design  purposes , 
to  ensure  that  excessive  momentum  pressure  drops  do  not  occur. 

It  should  be  emphasised  that  these  unusual  property  variations  and  the 
associated  transfer  phenomena,  are  confined  to  a  rather  narrow  range  of 
pressure  above  the  critical  pressure.  In  rocket  engines  where  combustion 
occurs  at  highly  supercritical  pressures  (P/P0  »  5)  these  conditions  will  not 
arise,  and  for  practical  purposes  the  fluid  can  be  considered,  above  a  certain 
temperature,  to  have  monotonically  varying  properties. 


2.  PROPERTY  VARIATIONS  IN  THE  NEAR-CRITICAL  REGION 

In  a  number  of  technological  applications,  severe  physical-property 
gradients  (both  radial  or  axial)  may  occur,  makinv  unreliable  the  use  of 
empirical  correlations  or  theoretical  analyses  based  on  the  constant  property 
conditions.  These  property  gradients  are  most  visually  associated  with 
gradients  of  temperature,  or  high  ratios  of  wall  temperature  to  fluid  bulk 
temperature,  since  all  the  relevant  physical  properties  of  fluids  are  to  some 
extent  temperature  dependent. 

/A  familiar  ..... 


A  familiar  example  of  this  is  the  modification  to  the  forced  convection 
correlation 


Mu  s  f  (fie,  Pr) 


1 


made  by  Sieder  and  Tate  (2)  to  account  for  the  effect  of  the  large  difference 
in  viscosity  between  the  bulk  fluid  and  that  at  the  wall  when  heating  or  cooling 
viscous  liquids,  finding  wide  application,  for  example,  in  commercial  tube  and 
shell  heat  exchangers: 


%  •  f  (®v  Sr)  . * 

In  flow  of  gases  at  supersonic  velocities,  extreme  wall-  to  bulk  fluid- 
temperature  ratios  are  encountered,  and  fluid  properties  vary  greatly  between 
the  cooled  wall  and  the  fluid  stream.  In  allowing  for  the  variation  of 
properties,  e.g.  across  the  fluid  stream  cross-section,  it  is  often  possible  to 
express  the  individual  transport  and  thermal  properties  as  simple  linear 
functions  of  temperature,  since  their  variation  is  monotcnic  with  temperature. 

A  similar  condition  may  occur  in  cases  where  heat  is  transferred  at  very  high 
heat  fluxes  to  liquids. 

The  region  in  which  pressure  is  greater  than  the  critical  pressure  of  the 
fluid,  while  the  temperature  may  be  dose  to,  but  below  or  above  the  critical 
temperature,  is  referred  to  here  as  the  "supercritical  region".  Of  special 
interest  is  the  region  of  pressure  above  the  critical  pressure,  up  to  reduced 
pressure  of  approximately  2,  where  small  changes  in  either  temperature  or 
pressure  can  lead  to  drastic  changes  in  the  magnitude  of  fluid  properties. 

This  region  is  referred  to  herein  as  the  "near-critical"  region. 

The  most  striking  variation  with  temperature  occurs  in  the  specific  heat 
at  constant  pressure,  which  at  the  critical  point  has  infinite  value.  Above 
the  critical  pressure,  the  C  isobars  each  have  a  pronounced  maximum  close  to 
the  critical  temperature.  Hie  temperature  at  which  this  maximum  value  occurs 
on  a  particular  supercritical  isobar,  referred  to  as  the  "transposed  critical 
temperature"  at  that  pressure',  increases  with  pressure,  the  value  of  the  maximum 
specific  heat  dedining  steadily  with  increasing  pressure. 

In  general  other  physical  properties  of  the  fluid  show  their  maximum  rate 
of  change  with  temperature  near  the  transposed  critical  temperature.  Measure¬ 
ments  of  the  thermal  conductivity  of  carbon  dioxide  (3)  and  ammonia  (4)  dose 
to  the  critical  point  hage  shown  that  marked  peaks  occur  in  this  region. 

Evidence  suggests  that  the  maxima  are  less  strong  and  attenuate  more  rapidly 
away  from  the  critical 
in  both  C  and  thermal 
density.  p 


point  than  those  of  the  specific  heat,  C  •  The  maxima 
conductivity  appear  to  occur  in  fluid  air  the  critical 
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2.1  Bolling-like  Phenomena  in  the  Near-critical  Region 

In  spite  of  these  drastic  property  variations  in  the  near-critical  region, 
the  fluid  is  usually  considered  to  he  a  monophase,  although  it  cannot  be 
considered  wholly  gaseous  or  liquid.  It  is,  however,  possible  to  consider  the 
existence  of  two  quasi-phases  above  the  critical  pressure,  under  transient, 
non-equilibrium  conditions,  which  could  occur  during  transfer  of  heat  to  the 
fluid.  By  an  examination  of  the  valuation  of,  for  example,  density,  viscosity, 
thermal  conductivity  and  specific  heat  of  the  fluid  in  the  near-critical  region, 
it  can  be  seen  that  the  transposed  critical  temperature  represents  a  point  of 
transition  from  the  type  of  property  characteristics  occurring  at  relatively 
low  temperatures,  to  that  typical  of  higher  temperatures. 

The  fluid  at  temperatures  well  below  the  transposed  critical  temperature 
exhibits  typical  liquid  property  characteristics:  the  specific  volume  variation 
indicates  a  low  value  of  thermal  expansion  coefficient}  viscosity  (ji )  and 
thermal  conductivity  (k)  are  high,  decreasing  markedly  with  temperature.  At 
temperatures  considerably  above  the  transposed  critical  temperature, 
conductivity,  viscosity  and  density  all  have  the  relatively  low  values  typical 
of  gases;  n  and  k  also  increase  slightly  with  temperature  and  the  fluid  density 
is  inversely  proportional  to  the  absolute  temperature. 

From  these  considerations,  Goldmann  (5)  postulated  a  mechanism  of  heat 
transfer  in  supercritical  fluids,  rather  similar  to  nucleate  boiling  and  based 
on  the  difference  in  physical  properties  between  the  dense  fluid  comprising  the 
bulk  of  the  fluid  stream,  and  the  "gas-like1*  fluid  close  to  the  wall  at  a  much 
higher  temperature .  There  is  experimental  evidence,  obtained  with  a  number  of 
fluids,  that  some  such  mechanism  occurs: 

(i)  peaks  in  the  local  heat  transfer  coefficient  distribution, 
related  to  surface  temperatures  (or  in  seme  cases  the  "film1*  temperature) 
close  to  the  transposed  critical  temperature  (6,7), 

(ii)  photographic  evidence  that  a  separation  of  pseudo-phases  occurs 
at  supercritical  pressures  in  pool  boiling  experiments  (8),  and 

(iii)  observations  of  severe  vibrations  or  oscillations  in  this  region, 
in  forced  convection  flow  systems  and  natural  circulation  loops,  for 
hydrogen  (9).  and  for  other  fluids  (p,l0,1l).  In  many  cases  these  effects 
have  been  recorded  simultaneously  with  large  increases  in  measured  heat 
transfer  rates  without  corresponding  increase  in  the  wall  temperatures,  a 
familiar  phenomenon  with  nucleate  boiling  in  flow  systems. 

Other  experimenters  have  not  observed  the  above  phenomena  associated  with 
"boiling-like"  behaviour.  Bather  the  properly  variations  of  the  near-critical 
region  have  been  manifested  in  unusual  axial  wall  temperature  distributions 
(12);  in  distinct  minima  in  local  heat  transfer  coefficients,  associated  with 
the  local  bulk  fluid  temperature  (13,1V,15)»  and  in  distorted  temperature  and 
velocity  profiles  across  the  flow  cress-sections  (l6).  In  a  number  of  these 
studies  in  the  critical  region  it  was  observed  that  under  certain  conditions, 
heat  transfer  coefficients  declined  consistently  with  increasing  heat  flux, 
e.g.  Beference  6,  as  is  predicted  by  the  theoretical  treatment  of  Deissler  (l7) - 
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The  occurrence  of  "bciling-Hke"  mechanism  is  evidently  dependent  on  the 
local  values  of  bulk  fluid  and  surface  temperatures,  and  their  relation  to  the 
transposed  critical  temperature. 

Il  general,  for  those  cases  where  surface  temperatures  of  the  same  order 
as  the  critical  temperature  were  attained,  but  the  bulk  fluid  was  much  below 
this  temperature,  phenomena  similar  to  nucleate  boiling  were  reported. 

In  the  other  group  of  experiments,  in  which  bulk  fluid  temperatures  close 
to  the  transposed  critical  temperature  were  encountered,  the  occurrence  of 
maxima  or  minima  in  the  local  values  of  heat  transfer  coefficients  wa3  related 
to  the  magnitude  of  the  temperature  difference  between  the  heated  surface  and 
the  fluid.  At  very  high  temperature  differences,  minimum  values  of  heat 
transfer  coefficient  occurred  when  the  bulk  temperature  was  close  to  the 
critical  temperature.  At  lower  temperature  differences,  the  local  heat  transfer 
coefficient  showed  a  maximum  value  under  these  conditions.  This  influence  of 
temperature  difference  suggests  a  further  analogy  with  boiling  phenomena. 

Nucleate  boiling,  occurring  at  low  temperature  differences,  sustains  high  rates 
of  heat  transfer,  but  at  much  higher  temperature  differences  is  superseded  by 
film  boiling  where  the  insulation  of  the  hot  surface  by  an  unbroken  vapour  film 
allows  only  much  reduced  rates  of  heat  transfer  to  the  bulk  of  the  liquid. 

There  is  strong  experimental  evidence  that  buoyancy  forces  become  important 
in  heat  transfer  at  near-critical  conditions,  and  may  account  for  a  number  of 
the  phenomena  described  above.  This  is  particularly  the  case  when  local  peaks 
in  surface  teaqperature  occur  at  locations  Mere  the  bulk  fluid  is  close  to  the 
transposed  critical  temperature*  Under  these  conditions  the  orientation  of  the 
tube,  and  also  the  direction  of  flow  in  vertical  tubes,  may  clearly  indicate 
the  influence  of  buoyancy  forces. 

Shitsman  (18)  investigated  the  occurrence  of  peaks  in  the  local  surface 
temperature  in  the  flow  of  near-critical  water  in  vertical  tubes,  and  showed 
that  when  the  flow  was  downwards,  i.e.  when  the  buoyancy  forces  were  opposed 
to  the  forces  associated  with  the  forced  flow,  the  peaks  did  not  occur.  This 
is  shown  in  Fig.  1,  which  is  reproduced  from  Reference  18;  this  also  shows 
that  the  surface  temperature  profiles  for  upward  flow  had  the  anomalous  peaks 
ab-ve  a  certain  heat  flux,  although  the  experimental  conditions  were  otherwise 
similar. 

Miropolskiy  et  all.  (19)  also  reported  the  occurrence  of  anomalous  peaks  in 
the  wall  temperature  for  upward  flow  of  near- critical  water,  as  shown  by  curve  3 
of  Fig.  2.  Curves  1  and  2  of  this  figure  also  indicate  that  the  surface  tempera¬ 
ture  profile  may  be  significantly  different  at  the  top  and  bottom  of  a 
horizontal  tubo. 

In  this  connection  may  be  mentioned  some  high-speed  photographic  studies 
made  in  this  Laboratory  of  heat  transfer  from  electrically  heated  rods  to 
keros  ine  flowing  in  a  rectangular  cross-section  channel.  These  confirm  that 
buoyancy  effects  are  also  important  for  the  case  where  T#  >  Tffi  »  T^.  These 
observations  were  made  at  a  reduced  pressure  of  approximately  1.5,  with  bulk 
temperatures  near  ambient,  and  surface  temperatures  somewhat  above  the 
estimated  critical  temperature  for  the  fluid  (398° C). 
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It  was  observed  that  with  turbulent  flow  of  near-critical  fluid,  the 
violent  activity  at  the  surface  was  predominantly  at  the  top  of  the  rod,  this 
effect  becoming  more  marked  as  the  flow  rate  was  reduced.  At  Reynolds  number 
of  approximately  2000  (taking  the  rod  diameter  as  characteristic  dimension) 
good  cine  films  were  obtained  at  up  to  4500  frames  per  second.  Very  violent 
activity  was  observed  to  occur  at  the  top  surface  only,  causing  turbulent 
density  changes  well  into  the  bulk  of  the  stream,  whereas  at  the  lower  surface 
there  was  virtually  no  activity.  This  can  be ^  seen  in  Fig.  3»  which  is 
reproduced  from  one  frame  of  the  relevant  cine  film. 

For  the  case  of  flows  through  vertical  tubes,  Brown  and  Cauvin  (20,21), 
in  studies  of  temperature  and  velocity  fluctuations  in  air  flowing  at  Reynolds 
numbers  between  270  and  6900,  have  demonstrated  that  turbulent  transport 
depends  on  the  interaction  of  free  and  forced  convection.  The  magnitude  of 
turbulent  heat  transfer  is  determined  by  the  temperature  and  radial  velocity 
fluctuations  in  the  fluid  stream,  and  the  correlation  coefficient  between  them* 
In  upward  flow  the  effect  of  buoyancy  is  to  reduce  the  correlation  coefficient, 
and  so  impair  the  effectiveness  of  the  radial  turbulent  fluctuations  in  the 
transfer  of  heat. 

Similarly,  in  near-critical  heat  transfer,  the  buoyancy  forces  in  forced 
upward  flow  can  produce  local  "laminarisation’1  of  the  flow,  which  is  manifested 
as  a  peak  in  the  wall  temperature. 
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erimental  Data:  Ranees  of  Variables  for  Supercritical  Crvoee 


Straight  Circular  Tubes 


Authors 


Ranges  of  Reduced  Variables 

tn _ •  a  f  o.. 


Pressure 


Szetela 

(50) 

miller  et  al. 

(25) 

Thompson  and  Geery 

(15) 

Core  et  al. 

(22) 

Perroud  et  al. 

(30) 

Hendricks  et  al. 

(33) 

Hendricks  et  al. 

(31) 

Wright  and  Walters 

(23) 

Hendricks,  Simone s.u 
and  Friedman 

(34) 

Aerojet 

(37) 

Bulk  Fluid 
Temperature 

<V 


Surface 

Temperature 

(T) 


Ratio 

TA 


Max 

He 

H 

(cal/a 


These  references  do  not  state  the  tube  orientation  explicitly,  t3 
given  here  being  deduced  from  figures  or  schematic  diagrams 


-  6  - 


TABLE  1 


Ranees  of  Variables  f  or  Supercritical  Cryogenic  Hydrogen 

Straight  Circular  Tubes 

Reduced  Variables 

Maximum 

Reynolds 

Number 

(Maximum) 

Tube 

Orientation 

k  fluid 
perature 

Surface 

Temperature 

<v 

Ratio 

TA 

Heat 

flux 

(cal/ cm2  sec) 

1.25 

to 

1.44 

3»6 

to 

13.2 

2.9 

to 

10.2 

69 

31  x  105 

Horizontal* 

0.86 

to 

3.1 

- 

- 

940 

130  x  10? 

Vertical 

Upward 

Flow* 

0.91 

to 

1.73 

30.5 

(Max) 

16.5 

(Max) 

313 

17.8  x  id5 

Vertical 

Upward 

Flow* 

0.93 

to 

2.5 

25.0 

(llax) 

6.3 

(Max) 

235 

- 

Vertical 

Upward 

Flow* 

D.76 

24.0 

(Max) 

- 

92 

Vertical 

Upward 

Flow 

5.9 

to 

5-3 

- 

94 

93  x  id5 

No 

information 

given 

5.66 

to 

2.0 

- 

15.0 

(Max) 

96 

- 

Vertical 

Upward 

Flow 

1.6 

to 

1-3 

- 

! 

34 

20  x  105 

Horizontal* 

1.0 

[to 

r3 

2.4 

to 

20.0 

11.0 

(Max) 

290 

40  x  105 

Vertical 

Upward 

Flow 

Lo  1  7.2 

f to  |  to 

1.6  ;  28.0 

8C0 

t 

- 

Vertical 

Upward 

Flow 

» 

f 


tot  state  the  tube  orientation  explicitly,  the  information 
teed  from  figures  or  schematic  diagrams' 
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The  ranges  of  experimental  data  relating  to  forced  convection  heat  transfer 
to  supercritical  hydrogen  flowing  in  straight  tubes  of  circular  cross-section 
have  been  summarised  in  Table  1  and  are  indicated  in  Fig*  4*  Data  for  curved 
and  non-circular  cross  section  tubes  are  given  in  Tables  3  (p*12)  and  3  (p.14). 


3*1  Straight  Circular  Tubes 

A  few  supercritical  data  were  obtained  by  Core  et  al*  (22)  in  the  region 
of  pressure  just  above  the  critical  (Pre^#  up  to  1*13)  end  with  mean  bulk  fluid 
temperatures  confined  to  a  narrow  range  close  to  the  critical  temperature 
(30  to  33. 5° K).  A  number  of  data  at  the  same  pressure  and  mass  flow  rate  in 
this  region  show  that  with  increasing  heat  flux,  the  local  heat  transfer 
coefficient  increases  to  a  maximum  value  and  declines  with  further  increase  in 
the  heat  flux.  The  range  of  bulk  fluid  temperatures  is  too  narrow  to  draw  any 
conclusions  as  to  the  reason  for  this  behaviour. 

Thompson  and  Geery  (15)  covered  wider  ranges  of  pressure  and  temperature, 
using  para-hydrogen,  although  reduced  pressures  below  3*6  were  not  studied. 

The  authors  concluded  from  their  experimental  data  that  heat  transfer  to  super¬ 
critical  cryogenic  hydrogen  occurred  in  two  distinct  regimes: 

Regime  "A" 

This  occurs  with  bulk  fluid  temperatures  between  the  critical  value  and 
about  80° K,  in  which  heat  flux  was  found  to  be  proportional  to  0.7  power  of  the 
temperature  difference  between  heated  surface  and  hulk  fluid. 

The  ''best-fit"  dimensionless  correlation  for  this  regime  was: 

®-0-64 
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There  was  wide  scatter  of  experimental  data,  but  the  majority  were  correlated  to 
±  305S  by  this  expression.  This  appears  to  be  an  unstable  region  of  heat 
transfer. 

Regime  "B" 

The  data  for  the  second  regime,  representing  more  stable  conditions,  with 
heat  flux  proportional  to  the  temperature  difference,  were  correlated  in 
similar  fashion  to  give: 


Nu 
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It  was  found  that  this  second  equation  cculd  also  be  used  to  correlate  data 
in  Regime  "A”  with  fair  accuracy  except  for  those  obtained  at  low  heat  fluxes 
and  low  T^T^- -ratios. 

/The . 


The  ranges  of  parameters  for  the  two  regimes  are  summarised  in  Tuble  2 


TABLE  2 

Data  of  Thompson  and  Geerv  (15)  ~  Ranges  of  Parameters 


REGIME  "A* 

REGIME  "B" 

Heat  Flux,  cal/ as?  sec 

4.3  to  192 

23.5  to  313 

Temperature  Difference,  surface 
to  fluid,  degK 

2.8  to  387 

117  to  809 

Heat  Transfer  coefficient, 
cal/ css?  sec  degK 

0.08  to  0.81 

Om  10  "to  0*47 

Temperature  Ratio,  T 

1.06  to  8.72 

3.45  to  16.5 

Reynolds  Number 

2.62  to  13.0 
x  10? 

2.96  to  17.5 

x  10? 

Prandtl  Number 

0.73  to  1.70 

0.62  to  1 .54 

Frtn  this  it  can  be  seen  that  there  is  no  satisfactory  criterion  to  determine 
which  mode  of  heat  transfer  will  occur  under  any  given  set  of  conditions.  The 
first  expression  above  agrees  quite  closely  with  the  equation  proposed  by 
Wright  and  Walters  (23)  to  correlate  their  data  for  gaseous  normal  hydrogen  at 
low  heat  fluxes;  while  the  relation  for  Thompson  and  Geery's  Regime  "B”  is  in 
agreement  with  the  equation  of  McCarthy  and  Wolf  (24)  for  heat  transfer  to 
gaseous  hydrogen  at  high  pressures  and  heat  fluxes. 


The  authors  *  claim  that  their  data  show  the  same  minimum  in  the 
generalised  heat  transfer  coefficient  (=-———)  as  reported  by  Powell  (l  3)  for 
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cryogenic  oxygen  at  reduced  pressure  of  1*4  does  not  seem  to  be  justified.  The 
minima  in  these  experiments  with  hydrogen,  at  reduced  pressures  between  3*6  and 
7*2,  are  associated  with  the  transition  from  their  Regime  "A"  to  Regime  "B", 
rather  than  with  the  dependence  of  fluid  properties  on  temperature  near  the 
critical  tesperature.  The  phenomenon  described  is  virtually  the  same  at  reduced 
pressures  cf  3*7  end  7*  At  these  pressures  the  tesperature  dependence  of  fluid 
properties  is  relatively  slight,  even  in  the  vicinity  of  the  critical 
tesperature. 


Miller  et  ml.  (23,26)  give  experimental  data  for  turbulent  forced 
convection  for  cryogenic  hydrogen  flowing  in  vertical  heated  tubes,  over  a  wide 
range  of  conditions.  Pressures  were  between  27  and  170  atm,  Reynolds  numbers 
(bulk)  up  to  13  x  10s  and  the  maximum  heat  flux  was  940  cal/ as?  sec.  The  data 
were  correlated  by  the  equation  of  Hess  and  Runs  (q.v.)  which  includes  a 
correction  for  the  wall-to-buik  fluid  kinematic  viscosity  ratio. 


-  8  - 
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At  the  lower  supercritical  pressures  studied  (reduced  pressure  of 
approximately  2),  large  increases  in  tube  wall  temperatures  were  observed  at 
locations  where  the  bulk  fluid  temperature  was  slightly  above  the  critical 
temperature  ("film-boiling"  type  of  effect  giving  a  minimum  in  the  local  heat 
transfer  coefficient).  This  effect  was  not  observed  at  higher  pressures  when 
cryogenic  hydrogen  behaved  in  a  similar  way  to  hydrogen  gas  at  higher 
temperatures. 

In  France,  Perroud  and  co-workers  have  published  a  number  of  papers 
(27,28,29,30),  giving  experimental  data  for  forced  convection  heat  transfer  to 
liquid  hydrogen  in  circular  tubes,  and  also  in  rectangular  tubes  in  which  only 
one  side  was  heated.  The  number  of  data  in  the  supercritical  region  is 
restricted  to  the  pressure  range  13  to  16  atm  and  for  circular  tubes  only. 

Fluid  (6C$i  .para-hydrogen)  entered  the  tubes  at  25°  K,  and  surface  temperatures 
up  to  800° K  were  attained.  In  this  region  of  pressure  and  surface  temperature, 
sharp  increases  in  wall  temperature  were  commonly  encountered  at  points  where 
the  local  fluid  temperature  was  in  the  close  vicinity  of  the  critical 
temperature.  Local  Nusselt  numbers  for  this  region  are  ouch  below  those 
predicted  by  the  correlations  recommended  in  the  above  papers  for  slab  critical 
conditions,  in  which  the  physical  properties  of  the  fluid  are  evaluated  at  the 
"film"  temperature,  viz.  T^  =  4r  (TJ}  T^).  Sapid  fluctuations  in  local  wall 
temperature,  pressure,  and  mass  flew  rate  were  also  observed  under  these  super¬ 
critical  conditions,  but  attenuated  on  further  increase  in  the  fluid  pressure. 

Experimental  data  on  heat  transfer  to  near-critical  cryogenic  hydrogen  are 
included  in  the  NASA  technical  note  of  Hendricks,  Graham,  Hsu  and  Friedmann 
(31),  in  which  they  are  compared  with  data  obtained  for  flow  of  gaseous  and 
two-phase  fluid.  Vertically  positioned  circular  tubes  were  used,  of  between 
0.476  cm  and  1.287  cm  internal  diameter,  and  length  40  to  60  cm.  The  fluid 
entered  the  tubes  at  temperatures  between  22°  K  and  67° K  (only  the  upward  flow 
direction  was  investigated)  and  the  range  of  pressures  covered  was  5*5  to  55 
bars. 


A  similarity  was  noted  in  axial  wall  temperature  distributions  for 
sub  critical  fluid,  and  for  supercritical  conditions  where  the  fluid  bulk 
temperature  at  the  inlet  was  much  below  the  transposed  critical  temperature 
(Tb  «  Tg).  In  both  these  cases  the  wall  temperature  first  rose  to  a  maximum 
and  then  declined.  There  was  also  a  similarity  between  the  observed  pressure 
drops  and  in  the  occurrence  of  fluid  oscillations  for  the  two  fluid  states. 
Fluid  oscillations  occurred  over  limited  "unstable"  regions  of  heat  flux  and 
fluid  flow  rate. 

The  axial  position  of  the  peak  in  surface  temperatures  occurred,  as  a 
crude  approximation,  where  the  bulk  temperature  approached  the  transposed 
critical  temperature,  but  also  depended  sensitively  on  heat  flux,  flow  rate 
and  the  enthalpy  of  the  fluid  at  the  inlet  to  the  heated  tube.  No  functional 
dependence  of  surface  temperature  on  these  or  other  parameters  could  be 
demonstrated. 

Measured  heat  transfer  coefficients  were  compared  with  valui  s  predicted 
by  the  Nusselt  film  correlation  for  gaseous  hydrogen  (as  given  in  kef.  32), 
when  it  was  found  that  the  ratio  of  these  two  coefficients  approached  unity 
as  the  pressure  increased.  For  example,  in  a  tube  of  diameter  0.89  cm  and  at 
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axial  location  11.25  cm  from  the  inlet  fiange(ir  =  12.6)  the  ratio  of  experi¬ 
mental  to  calculated  film  coefficient  was  found  to  have  values  between  3  ant'  4 
at  reduced  pressure  1.2  (i.e.  13*8  bars);  at  55  bars,  or  reduced  pressure  4*2, 
the  ratio  fell  to  approximately  1«5» 

This  reference  also  includes  experimental  data  for  the  condition  of  axial 
heat  flux  gradient,  which  was  produced  by  using  an  electrically  heated  tube 
with  tapering  wall  thickness,  giving  the  maxi  muni  power  generation  at  the  middle 
of  the  tribe.  It  was  found  that  the  minimum  heat  transfer  coefficient  then 
coincided  with  the  position  of  maximum  heat  flux;  also  that  the  relative 
reduction  in  heat  transfer  coefficient  here  was  considerably  more  marked  in  the 
case  where  inlet  bulk  fluid  temperature  was  low  (2?.3°K  against  37«1°X).  In  the 
case  of  low  inlet  bulk  temperature  the  initial  value  of  heat  transfer  coefficient 
was  considerably  higher,  and  in  this  case  the  bulk  fluid  temperature  approached 
the  transposed  critical  value  in  the  high-heat  flux  region  of  the  tube. 

Various  methods  of  correlation  are  discussed  in  this  note  in  terms  of  the 
Martinelli  two-phase  parameter,  which  it  had  been  demonstrated  earlier  (33) 
could  be  reformulated  and  used  to  correlate  near-critical  heat  transfer  data. 

The  essential  features  of  this  method  of  correlation  are  discussed  in  Section  4 
below  (page  22).  In  this  more  recent  note  the  method  is  considerably  modified 
by  inclusion  of  the  Sterman  parameter  ( Sr)* for  supercritical  fluid*  anddann 
entrance  effect  term;  the  experimental  Nusselt  number  is  based  on  enthalpy 
differences. 

The  new  correlation  involves  considerable  computational  effort  and  there 
is  gross  scatter  of  experimental  data  over  the  large  range  of  supercritical 
pressures  covered  (13*8  -  55  bars).  In  common  with  many  other  proposed  super¬ 
critical  hydrogen  heat  transfer  correlations,  the  amount  of  scatter  can  be 
reduced  to  acceptable  proportions  (in  this  case  within  +35  to  -25$  if  data  in 
the  extended  entrance  section  of  the  tube  (i/d  <  JO)  are  excluded. 

The  range  of  high  supercritical  pressures  (69  -  172.5  bars)  wa3  covered 
in  experiments  reported  in  another  NASA  report  by  Hendricks,  Simone au  and 
Friedman  (34)*  Liquid  hydrogen  entered  at  between  33  and  110°K  and  flowed 
upwards  through  vertical,  uniformly  heated  tubes  with  diameters  between  0.536 
and  1.114  un.'  Surface  temperatures  were  between  00  and  670° K,  and  the  maximum 
heat  flux  was  290  cal/cm2  sec.  Behaviour  of  the  fluid  in  this  range  of  pressure 
was  completely  gas-like,  and  the  majority  of  the  data  could  be  correlated 

/using . 


*The  Sterman  number  (Sr),  originally  proposed  for  use  in  flow  with  boiling, 
through  heated  tubes,  is  the  ratio  of  the  transverse  to  axial  velocities,  and 
a  measure  of  the  penetration  of  the  bulk  fluid  stream  by  the  vapour  generated 

at  the  hot  wall.  In  two-phase  flow,  jS-  may  be  regarded  as  the  mean  transverse 
velocity  relative  to  the  total  wall  area. 

For  supercritical  conditions,  the  Sterman  number  becomes  where 

b 

A7  and  AH  are  respectively  the  enthalpy  and  specific  volume  differences 
between  the  hypothetical  light  ("ges-like")  and  heavy  ("liquid-like")  species 
referred  to  on  page  22. 


using  tha  Hus salt  fliui  type  ox  relation:  Nu^,  =  u  ne^*  Pr^° ' 4 ,  where  the 
exponent  p  has  a  value  between  0.8  and  0.9,  depending  on  tube  diameter. 

Oscillations  accompanied  by  a  significant  increase  in  heat  transfer 
coefficient  were  observed  under  certain  restricted  conditions  with  the  largest 
diameter  tubes  employed  (1.114  cm)  only.  The  effect  of  pressure  on  this 
phenomenon  is  not  discussed  but  oscillations  appear  to  have  been  confined  to 
the  range  69  -  103.5  bars,  and  there  is  evidence  that  they  were  affected  by 
heat  flux.  However  similar  acoustic  vibrations  were  also  observed  in  tht-< 
apparatus  when  there  was  no  heating,  and  this  always  occurred  at  very  low  flew 
velocities. 

3.2  Curvilinear  Tubes  and  Non-circular  Cross  Sections 

The  following  work  includes  experimental  data  obtained  with  tubes  other 
than  straight  ones  with  circular  cross-section  and  are  summarised  in  Table  3* 

Hendricks  and  Simon  (35)  investigated  heat  transfer  to  hydrogen  including 
supercritical  pressures  (20.7  and  41 *4  bars)  in  curved  tubes  of  various 
diameters  and  geometries.  The  tubes  were  vertical,  and  uniformly  heated  by 
electric  power;  the  apparatus  was  a  "blow-down"  loop.  Many  of  the  data  were 
apparently  obtained  at  10.35  bars  (reduced  pressure  0.8)  and  are  referred  to 
in  discussion  as  "near-critical":  it  is  with  these  data  that  the  most  pronounced 
effects  of  tube  curvature  are  reported.  Heat  transfer  was  improved  at  concave 
surfaces  (viewed  from  within  the  tube)  and  diminished  at  convex  surfaces.  The 
ratio  of  concave  to  convex  heat  transfer  coefficients  exceeded  3  for  data  at 
reduced  pressure  0.8,  but  the  magnitude  of  the  ratio  depended  on  the  radius  of 
curvature ,  angular  position,  tube  diameter  and  fluid  structure.  The  influence 
of  tube  diameter  makes  generalisation  of  these  effects  impossible,  since  heat 
transfer  coefficients  at  both  the  convex  and  concave  surfaces  may  either 
maximise  or  minimise  depending  up^n  the  tube  diameter.  The  ratio  kgongAcouy 

has  a  minimum  value  (for  1.77  cm  diameter  curved  tubes)  near  the  point  in  the 
tube  where  fluid  bulk  temperature  equals  the  transposed  critical  temperature; 
for  0.8  cm  diameter  tubes  with  11.1  cm  radius  of  curvature,  however,  the  ratio 
remained  constant. 

Contrary  to  the  authors'  statement  the  graphical  data  for  supercritical 

hydrogen  presented  indicate  that  at  this  pressure  (41  *4  bars)  heat  transfer 

coefficients  at  both  concave  and  convex  surfaces  are  slightly  higher  than  those 

predicted  by  the  Dittus-Boelter  type  of  equation  for  straight  tubes.  From  this 

graph  the  ratio  h  Jh _ has  values  between  1  and  1.2  over  the  length  of  a 

cone  conv 

0.95  cm  diameter  tube  with  radius  of  curvature  19  cm. 

It  seems  likely  that  both  the  centrifuging  and  secondary  flow  effects 
reported  for  reduced  pressure  0.8  may  also  occur  in  curved  tubes  at  pressures 
immediately  above  the  critical  pressure*  (it  will  be  noted  that  supercritical 
pressures  below  reduced  pressure  1.6  are  not  covered  in  the  present  reference.) 
However  there  is  some  ambiguity  in  the  authors'  discussion  of  their  results, 
because  of  their  omission  to  make  clear  the  pressure  conditions  for  the  data 
presented  and  discussed. 
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Thompson  and  Ceery  (36)  carried  out  an  analysis  of  heat  transfer  dat6  for 
liquid  hydrogen  flowing  in  both  curved  and  straight  tubes.  [This  report  is 
unobtainable  in  the  UK  at  present  and  details  are  quoted  from  the  later  Aerojet 
report  (37)  of  1%7»  It  is  not  clear  from  this  whethor  any  original  dp-*'*  are 
given  in  Reference  36.]  They  concluded  that  the  data  for  pressure 
41*4  and  103*5  bars  were  best  represented  by  the  following  relati«- 

v 

Nuf  =  0.0208  Ref°*8  Prf°'4  (l  +  0.01457~)  ■  •  <3A 

'b 

This  is  the  expression  of  Hess  and  Runs  (30)  discussed  below  (p.  21 ),  with 
inclusion  of  CT ,  which  is  a  function  of  bulk  fluid  temperature,  as  follows: 

If 


(The  influence  of  the  radius  of  curvature  on  the  data  is  not  given  in  Ref.  37) 


Since  the  enhancement  of  CL  at  low  coolant  temperatures  is  apparently  not 
a  function  of  pressure,  this  influence  of  fluid  temperature  cannot  be  the 
result  of  fluid  property  variations  in  the  region  of  transposed  critical 
temperature,  (in  the  range  of  reduced  pressure  covered,  3*2  to  8.0,  these  will 
in  any  case  be  relatively  slight.)  It  will  be  noted  that  the  correlation 
suggested  indicates  that  the  ratio  of  heat  transfer  coefficients  for  curved  to 
straight  tubes  is  constant  over  the  range  of  bulk  fluid  temperature  quoted. 

A  second  Aerojet  REON  report,  also  not  available  in  UK  at  the  time  of 
writing,  and  dating  from  Hay  1963  \3S)  I  is  also  mentibned'T>rie£ly  in  Ref.  37* 
Heat  transfer  to  hydrogen  at  near  critical  temperatures  and  supercritical 
pressures,  and  flowing  turfculently  in  straight  and  curved  tubes,  was  studied 
in  a  one-tube  geometry.  For  curved  tubes  with  non-uniform  heating  some 
improvement  in  heat  transfer  is  reported,  but  it  was  not  possible  to  generalise 
the  results  in  a  form  permitting  extrapolation  to  other  design  conditions. 

/The . 
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The  1%7  NASA  report  (37)  ^  giving  the  results  of  Aerojet  studies  in 
support  of  the  Phoebus  nozzle  development  programme,  provides  the  most  compre¬ 
hensive  collection  of  data  fvom  US  sources  on  heat  transfer  to  supercritical 
hydrogen,  under  conditions  of  non-uniform  heating  and  for  non-circular 
geometries.  The  apparatus  is  a  Mbiow-downM  loop,  described  in  Reference  39* 

The  study  includes  data  for  straight  tubes  of  length  7»5  and  15  cm  and 
diameter  0.373  cm,  and  uniform  nail  thickness.  These  tubes  were  used  to 
determine  the  effects  of  coolant  temperatures  close  to  the  transposed  critical 
temperature  on  local  heat  transfer  coefficient,  and  also  to  demonstrate  that 
high  heat  fluxes  of  the  order  580  cal/cm2  sec  can  be  sustained  when  using 
supercritical  hydrogen  as  coolant. 

The  straight  tube  data  are  compared  with  the  equation  of  Hess  and  Kiunz 
(see  p#  21).  This  underestimates  the  experimental  data,  which  show  considerable 
scatter,  by  up  to  7Q$  from  the  correlating  equation,  but  the  inclusion  of  u 
kinematic  viscosity  ratio  in  the  relation  accommodates  widely  differing  wall 
temperature  regimes.  It  is  attempted  to  separate  the  entrance  effects  from 
bulk  temperature  effects  in  the  short  tubes.  Both  trends  are  in  fact  slight 
compared  with  che  amount  of  scatter  of  experimental  data,  although  heat  transfer 
coefficients  tend  to  be  higher  at  low  bulk  temperatures. 

The  ranges  of  tube  geometry  parameters,  and  of  operating  conditions  for  the 
experiments  with  curved  tubes  are  summarised  in  Table  5* 

TABIB  5 

Summary  of  Ranges  of  Test  Parameters  for  Curved -Tub  63  (Ref.  37) 


Parameter 

Circular 

cross-section  tubes 

Non-circular 
cros3-section  tubes 

R||  cm 

8.2 

4.2 

11.3 

5'0 

R^  cm 

16.4 

8.2 

18.1 

8.2 

Dg  cm 

0.44 

0.44 

0.40 

O.40 

I/De 

(heated  length) 

21-69 

23  -  47 

23  -  69 

26  -  51 

P  bars 

59  -  87 

61  -  81 

75  -  87 

60  -  85 

\  °K 

36-52 

40  -  62 

44-72 

39-6 

T  °K 

3 

108  -  912 

145  -  890 

226  -  749 

93  -  595 

&  g/ cm2  sec 

579  -  894 

388  -  748 

420  -  477 

479  -  825 

q  cal/c®2  sec 

215  -  361 

186  -  375 

210  -  364 

203  -  375 

/Figure  1 
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Figure  5  indicates  the  form  of  the  tubes  and  their  cross  sections.  The  curved 
tubes  simulated  two  different  nozzle  throat  geometries,  and  were  used  to 
determine  the  degree  of  enhancement  of  heat  transfer  coefficients  caused  by 
tube  curvature.  Ncn-asymmetric  neat  fluxes  were  obtained  by  variation  of  wall 
thickness  about  the  periphery  of  both  circular  and  non-circular  tubes  (the 
latter  were  flattened  circular  tubes  with  uniform  flow  area,  as  shown  in  cross 
section  in  Figure  5)« 

The  curved  tubes  had  long  entrance  lengths  of  over  20  diameters  to  avoid 
any  entrance  or  coolant  temperature  effects.  This  allowed  the  establishment  of 
the  thermal  boundary  layers,  and  preheating  of  the  coolant  to  a  temperature  of 
approximately  45°  K  before  the  fluid  reached  the  curved  portion  of  the  tubes.  A 
comparison  of  the  experimental  data  with  the  Hess  and  Kunz  relation  reveals  the 
folio. ring  points: 

(a)  heat  transfer  coefficients  at  the  concave  surfaces  (i.e.  outside  of 
the  curve)  were  enhanced  by  factors  up  to  2:1, 

(b)  the  magnitude  of  this  enhancement  decreased  when  the  radius  of 
curvature  was  increased, 

(c)  neither  asymmetrical  heating  nor  non-circular  ducts  caused  the 
measured  local  heat  transfer  coefficients  to  be  lower  than  predicted 
by  the  Hess  and  Kunz  relation,  and 

(d)  at  those  positions  of  the  curved  tube  surfaces  corresponding  to  the 
highest  heat  flux  encountered  in  a  rocket  nozzle,  the  measured  heat 
transfer  coefficients  were  between  1.0  and  1.4  times  the  values 
predicted  by  the  correlating  equation  of  Hess  and  Kunz. 

The  significance  of  the  results  with  curved  tubes  for  the  design  of  the 
Phoebus-2  nozzle  is  discussed.  Friction  factors  were  investigated  and  tend  to 
confirm  the  heat  transfer  results,  though  the  findings  are  not  conclusive. 

The  occurrence  of  vibrations  was  noticed  in  some  tests,  but  it  was  not  possible 
to  demonstrate  any  effect  on  heat  transfer. 

3*3  Other  Geometries 

A  number  of  studies  have  provided  data  on  heat  transfer  to  hydrogen  in 
flow  ducts  other  than  electrically  heated  tubes  in  test  rigs,  such  as  heat 
exchangers,  and  these  provide  additional  information  on  phenomena  peculiar  to 
the  supercritical  state. 

A  pertinent  reference  to  some  observations  made  under  conditions  of  pool 
heating  is  also  included  in  this  Section. 

Thurston  (9)  studied  the  occurrence  of  pressure  oscillations  during  heat 
transfer  to  liquid  hydrogen  in  an  apparatus  which  cooled  from  ambient  to 
cryogenic  temperatures  during  each  run.  The  nature  of  the  acoustic  modes  of 
vibration  observed  depended  very  much  on  the  state  of  the  fluid.  In  the  super¬ 
critical  condition  both  Helmholtz  and  open-pipe  resonance  were  noted,  also  & 
"supercritical"  mode  which  occurred  only  when  pressures  exceeding  the  critical 
value  werer.Atfcained.  None  of  these  effects  was  observed  when  the  liquid 

/hydrogen  . 
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hydrogen  entered  the  test  section  at  temperatures  abcve  the  transposed  critical 
temperature,  from  which  it  is  concluded  that  the  presence  of  dense  "liquid- 
like"  fluid  is  a  necessary  condition  for  these  vibration  effects  to  take  place 
during  heat  transfer  to  supercritical  hydrogen. 

Bartlit  and  V/illicmson  (40,41 )»  studied  transfer  of  heat  between  r/ater  at 
near-ambient  temperatures  and  cryogenic  hydrogen  in  simple,  one-pass  annular 
heat  exchangers.  It  was  found  that  overall  heat  transfer  coefficients  as  low 
as  l/2  to  l/9  of  those  predicted  by  the  dimensionless  correlation  of  Diitus 
and  Boelter  (42)  were  obtained,  unless  a  swirl-inducing  strip  was  placed  in  the 
upstream  end  of  the  exchangers,  when  results  in  reasonable  agreement  with  the 
Dittus  and  Boelter  relation  were  obtained. 

These  results  are  explained  by  a  separation  between  a  core  of  dense 
hydrogen  flowing  along  the  centre  of  the  tube,  and  the  warmer  fluid  layer  at 
the  walls.  The  effect  of  pressure  is  not  discussed  by  the  authors,  and  it  is 
evident  from  the  tabulated  data  that  this  phenomenon  is  not  a  consequence  of 
supercritical  property  variations,  for  all  the  experiments  in  which  it  was 
observed  were  conducted  at  subcritical  pressures.  The  phenomenon  described  is 
consistent  with  the  occurrence  of  film  boiling,  which  would  be  expected  at  the 
large  overall  values  of  temperatures  difference  between  the  water  and  the 
cryogenic  fluid,  with  hydrogen  at  subcritical  pressures.  It  is  not  clear 
whether  experiments  were  also  conducted  with  hydrogen  at  supercritical  pressures 
and  without  a  swirl  inducer  in  the  apparatus;  nor  whether  core-separation 
effects  were  observed  under  these  conditions.  It  was  noted,  however,  that  for 
all  runs  carried  out  without  a  swirl  inducer  in  the  exchanger,  large  pressure 
oscillations  were  observed. 

Graham  et  al.  (43)  studied  the  effects  of  acceleration  on  pool  heating  of 
liquid  hydrogen  at  both  subcritical  and  supercritical  pressures.  From  the 
evidence  of  the  high  speed  photographs  obtained  and  from  the  experimental  data 
the  authors  concluded  that  the  mechanisms  of  heat  transfer  in  established  film 
boiling  and  in  supercritical  heating  (respectively)  are  very  similar.  At  high 
temperature  differences  between  surface  and  fluid*  the  data  for  supercritical 
conditions  and  for  film  boiling  tended  to  merge  irrco  one  bend;  and  ell  super¬ 
critical  data  were  characterised  by  low  values  of  the  heat  transfer  coefficient. 

The  effect  of  increased  gravity  (7  g)  on  supercritical  heating  of  para- 
fcydregen  was  to  increase  rates  of  heat  transfer  by  approximately  65/2. 


4*  ANALYSIS  OF  SUPERCRITICAL  H5AT  TRANSFER  AND  CORRELATION  OF  EXFERIKMTAL  DATA 

A  number  of  theoretical  and  semi-theoretical  treatments  of  the  problem  of 
convective  heat  transfer  with  variable  physical  properties  has  been  published. 

In  some  cases  correlating  equations  have  resulted;  a  number  of  these  are 
intended  for  design  purposes  with  liquid  hydrogen  systems,  and  experimental 
data  were  employed  in  their  derivation. 

Deissler's  analysis  was  initially  developed  (44)  to  deal  with  the  case  of 
heat  transfer  to  turbulent  flow  of  gases  with  very  high  heat  fluxes,  where  the 
Prandtl  number  was  near  unity  and  constant  specific  heat,  C  ,  corid  be  assumed 
without  serious  error.  The  case  was  further  simplified  in  ?hat  viscosity  and 
thermal  conductivity  could  be  represented  as  comparatively  simple  functions  of 
temperature. 


The  analysis  is  based  on  the  analogy  between  transfer  of  heat  and  momentum 
exchange  in  turbulent  flow  and  employs  familiar  equations  for  molecular 
conduction  of  heat,  and  transfer  of  momentum  by  shear  stress,  modified  to  take 
account  of  the  additional  transfer  due  to  the  turbulent  motion  of  eddies  in 
the  fluid,  viz. 
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It  was  assumed  that  the  eddy  di^fusivities  for  heat  and  momentum  were 
equal:  =  ea;  and  further  that  the  expression  developed  for  eddy  diffusivity 

for  adiabatic  flow  conditions  can  be  applied  to  transfer  of  heat  with  variable 
fluid  properties. 

Deissler  employed  two  expressions  for  the  eddy  diffusivity  at  a  point  in 
the  fluid  stream;  one  for  locations  dose  to  the  wall,  viz. 


the  second  for  points  at  a  distance  from  the  wall,  vis. 
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in  which  both  n  and  m  were  assumed  to  have  constant  values,  derived  from 
experimental  data  for  adiabatic  flow  of  air. 

Velocity,  temperature  and  distance  from  the  wall  were  represented  as 
dimensionless  variables,  with  the  dimensionless  heat  flux  as  an  independent 
parameter.  Velocity  and  temperature  profiles  across  the  fluid  stream  were 
calculated  by  solution  of  the  resulting  simultaneous  differential  equations 
and  these  in  turn  were  used  to  calculate  heat  transfer  coefficients.  The 
method  involved  tedious  calculations,  and  from  a  number  of  such  calculations, 
Deissler  and  Eian  (45)  were  able  to  correlate  experimental  results  for  sir  to 
give  Nusselt  number  as  a  single  function  of  the  Reynolds  number,  when  fluid 
properties  were  evaluated  at  a  reference  temperature,  T  ,  where: 


*  \  +  1  -  V 
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with  x  =  0.4. 


Deissler  later  developed  this  analysis  to  cover  the  case  of  heat  transfer 
to  turbulent  flow  of  supercritical  water  (17),  which  is  considerably  more 
complex  since  Pr  >  1,  and  because  of  the  unusual  property  variations  in  the 
critical  region,  with  strong  peaks  in  the  specific  heat,  C  ,  and  in  the 
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Fr-ttfidtl  number,  Pr ,  close  to  the  critical  tenperature  of  the  fluid.  A3 
previously,  profiles  of  temperature  and  velocity  across  the  fluid  stream  were 
obtained,  but  because  of  the  nature  of  fluid  property  variations  with 
temperature ,  these  profiles  had  to  be  calculated  separately  for  each  wall 
temperature  initially  assumed.  The  results  show  a  complex  dependence  of 
Nusselt  number  on  the  dimensionless  heat  flux  parameter,  /3.  In  order  to 
eliminate  the  effect  of  /?  it  was  necessary  to  determine  fluid  properties  at  a 
reference  temperature  which  was  a  function  of  both  the  wall  temperature  and  the 
ratio  of  wall  to  bulk  fluid  iemoiratures 

In  his  discussion  of  Deissler's  paper  (l7)»  Eckert  shows  a  simple  method 
of  overcoming  this  difficulty  by  plotting  the  reference  temperature  as  a 
function  of  the  transposed  critical  temperature,  T^.  This  effectively  reduces 
the  data  for  different  wall  temperatures,  T  ,  to  a  single  curve,  from  which  it 
is  seen  that  the  value  of  the  reference  temperature  to  be  used  is  largely 
determined  by  the  transposed  critical  temperafcire,  T^.  It  is  also  noted  that 
in  cases  when  the  bulk  fluid  temperature  is  close  to  T  the  correct  reference 
temperature  i3  T^;  similarly  that  the  surface  temperature,  T^»  should  be  used 
as  the  reference  temperature  when  T3  is  close  to  T  •  Table  6"  gives  the 
reference  temperature ,  T  ,  to  be  used  for  different  values  of  the  function  of 

V  ^s  and  V 

TABLE  6 


T 

X 

<  0 

*b 

0  to  1.0 

Tm 

>  0 

T 

s 

Goldmann  (46}  used  a  similar  analytical  method  to  Deissler,  but  redefined 
the  dimensionless  wall-distance  and  velocity  parameters  to  take  account  of  the 
assumption  that  the  eddy  diffusivity  at  a  point  in  the  flow  field  is  a  function 
of  the  fluid  properties  at  that  point,  but  is  not  affected  by  small  variations 
close  to  the  point.  The  method  allows  the  calculation  of  enthalpy  and  velocity 
profiles  from  an  assumed  set  of  wall  conditions,  vis.  T  ,  qs  and  r  ,  and 
correlation  of  experimental  data  is  in  terms  of  a  generalised  heat  flux 
parameter: 
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-  Hail,  Jackson  and  Khan  (47)  used  Goldmann’s  basic  assumption  but 
extended  the  analysis  to  take  account  of  the  effect  of  the  greatly  increased 
magnitude  of  the  expansion  coefficient: 


i  (&\  = 

7  v  Vp 


=  f  (t) 


of  the  fluid  close  to  the  heated  wall  when  this  is  close  to  the  transposed 
critical  temperature.  They  supposed  that  this  must  enhance  the  turbulent 
diffusivity  in  the  wall  region  in  proportion  to  the  degree  of  thermal  expansion 
in  the  fluid,  which  sharply  maximises  at  the  transposed  critical  temperature. 

The  modified  theory  shows  reasonable  agreement  with  the  trends  of  their  experi¬ 
mental  data  which  were  obtained  with  turbulent  flow  of  supercritical  carbon 
dioxide  between  parallel  plates,  under  conditions  of  fully  developed  velocity' 
and  temperature  profiles  across  the  flow.  The  model  proposed  here  is  in 
essence  the  "boiling-like"  mechanism  of  Goldmann  (5),  but  use  of  the  maxima  in 
thermal  expansion  coefficient  allows  this  effect  to  be  expressed  quantitatively. 
It  is  interesting  to  note  that  agreement  between  experiment  and  theory  above 
is  improved  by  introducing  a  scaling-factor  to  the  turbulent  diffusivity 
distribution  in  the  sub-layer  region  (both  distributions  are  consistent  with 
accepted  velocity  distribution  data  for  the  region).  The  scaling-factor  used 
by  Hall  et  al.  fell  linearly  from  a  value  of  20  at  the  wall  to  unity  near  the 
edge  of  the  sub-layer,  and  thus  represents  a  considerable  distortion  of  the 
original  distribution.  Figure  6  indicates  a  hypothetical  distribution  of  eddy 
diffusivity: 


=  y  Cl  *  exp  (~y)3 


where  e  is  in  arbitrary  units,  and  y  =  ^  is  the  dimensionless  wall  distance, 

a  y* 

y*  'jeing  the  thickness  of  the  laminar  sub-layer.  The  figure  also  shows  this 
distrib  ltion  multiplied  by  the  scaling-factor  employed  by  Hall  et  al.  In 
reality  the  diffusivity  will  also  depend  on  velocity  gradient  and  properties 
of  the  fluid  (notably  the  kinematic  viscosity)  as  functions  of  the  wall  distance 


i 


Nevertheless,  it  is  seen  that  the  modification  in  distribution  £  of  Fig.  6 
represents  a  significant  increase  in  the  influence  of  turbulence  eddies  within 
the  sub-layer,  which  is  required  to  enable  the  thermal  expansion  mechanism 
proposed  to  fully  account  for  the  increjises  in  heat  transfer  coefficients 
observed  experimentally.  Under  conditions  remote  from  the  critical  state  it  is 
evidently  possible  to  use  a  number  of  different  distributions  of  eddy  diffusivity 
within  the  sub-layer  which  will  be  consistent  with  measured  velocity  distri¬ 
butions,  and  will  also  enable  fair  agreement  with  heat  transfer  experiments  to 
be  obtained.  Uhen  near-critical  conditions  exist  at  the  wall  it  is  suggested 
that  the  analysis  is  far  more  sensitive  to  the  e  distribution  assumed,  since 
this  is  in  effect  "amplified"  by  the  nature  of  fluid  property  variations  within 
the  layer. 
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The  experiments  cf  Kali  et  al.  (47)  did  not  of  course  involve  conditions 
under  which  deterioration  of  local  heat  transfer  is  observed,  such  as  has  been 
described  for  certain  conditions  in  normal  pipe  flow.  Neither  does  their 
theory  outlined  above  predict  such  behaviour  which  has  been  discussed  earlier 
in  terms  of  the  influence  of  buoyancy  forces  (p.  4/» 

Evidence  of  periodic  heat  transfer  fluctuations  in  the  laminar  sub-layer 
is  given  in  a  recent  experimental  study  by  Armistead  and  Keyes  (48)  of  turbulent 
flow  of  unpressurised  water,  using  hot  film  sensors.  This  demonstrates  the 
existence  of  moderately  high  levels  of  turbulence  in  the  laminar  sub-layer,  and 
gives  further  cause  to  speculate  whether  the  eddy  diffusivity  distributions 
employed  in  many  theoretical  models  for  supercritical  heat  transfer  omit  some 
important  turbulence  promoting  mechanism  occurring  in  the  laminar  region. 


Hsu  and  Snith  (49) >  modified  Deissler’s  treatment  of  the  problem  to  take 
account  of  the  effect  of  the  variation  of  fluid  density  across  the  stream  on 
ths  local  values  of  the  eddy  diffusivities: 


where 
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The  procedure  for  calculating  velocity  and  temperature  profiles,  and  hence 
heat  transfer  coefficients  and  Nusselt  numbers,  was  similar  to  that  in  the 
analyses  already  discussed. 

For  flow  in  vertical  tubes,  the  influence  of  large  fluid  density 
valuations  in  inducing  natural  convection  in  near-critical  heat  transfer  was 
considered.  The  analysis  predicts  that  Nusselt  numbers  should  be  increased 
under  conditions  where  the  transposed  critical  temperature  lies  between  the 
wall  and  bulk-  fluid  temperatures.  The  effect  should  only  be  significant  for 
low  Reynolds  numbers  and  relatively  high  Grashof  numbers,  since  the  gravitation 
forces  are  increased  in  proportion  to  the  ratio  Gr/Re  . 

Szetela  (50)  obtained  experimental  data  for  hydrogen  at  bulk  temperatures 
between  33*3  and  47°K,  wall  to  bulk  temperature  ratios  from  3  to  11  and  reduced 
pressures  of  1.13  and  1.68.  A  limited  number  of  data  were  selected  and 
calculations  made  to  obtain  predicted  values  of  the  heat  transfer  coefficients: 
from  the  equations  of  Deissler  (17)  and  ilsu  and  Smith  (49) »  at  the  sane 
conditions.  Heat  transfer  coefficients  calculated  from  the  unmodified  equations 
of  Deissler  all  fall  below  the  experimental  values  by  differences  ranging  up  to 
335S.  The  coefficients  calculated  from  Hsr.'s  modified  equations  were  from 
50  to  1%  above  the  experimentally  measured  data. 

These  analyses  have  therefore  proved  inadequate  to  deal  with  the  extreme 
variation  of  fluid  properties  with  temperature  in  the  critical  region  for 
hydrogen;  also,  experiments  indicated  large  variations  in  temperature  along 
lie  tube  wall,  whereas  the  theoretical  treatments  discussed  above  assume 
uniform  wall  temperature.  It  also  proved  impossible  to  correlate  the  data 
using  the  conventional  dimensionless  relation  for  Nusselt  number,  with  fluid 
properties  evaluated  at  the  "film"  temperature,  T^.,  because  of  excessive 
scatter  of  the  data. 
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A  sssd.-65ipixu.eal  suuay  was  made  by  Hess  and  Kunz  (38) »  and  the  resoles 
compared  with  data  from  the  literature  for  hydrogen  in  the  supercritical 
region.  This  analysis  starts  from  the  same  basic  equations  as  those  previously 
discussed,  but  employs  a  single  expression  for  the  momentum  eddy  diffusivity: 


-rn  =  y2  [1  -  exp  rj))2  0 


the  factor  in  brackets  accounting  for  viscous  danping  of  turbulent  eddies  in  the 
region  close  to  the  wall.  A  comparison  was  made  of  the  computed  Nusselt 
numbers  with  experimental  values  derived  from  fieferences  33  and  50  and 
additional  unpublished  data.  At  low  heat  fluxes  good  agreement  was  obtained, 
but  with  increasing  heat  flux,  the  ratio  Ku  (expt)/Nu  (calc)  inoreased 
markedly.  It  was  found  that  the  increase  in  this  ratio  could  be  correlated  as 
a  linear  function  of  the  kinematic  viscosity  ratio,  'v^/v^  [v  »  h/p)»  W»a 
maximum  value  of  this  ratio  in  the  experimental  data  considered  was  140,  and  it 
is  evident  that  the  analysis  was  unable  to  account  for  the  variation  in  fluid 
properties  across  the  stream  under  such  extreme  conditions. 

Since  ti.*s  largest  variations  in  fluid  properties  occur  olose  to  the 
wall,  the  value  of  the  "damping  factor"  A4  in  the  non-dimension*!  form  of  the 
expression  for  eddy  diffusivity,  vis: 


jZ7pJ  =  f1  -  ex?  37  . ii 

was  examined  more  closely.  Values  of  A4  were  calculated  to  give  agreement 
between  the  experimental  and  computed  values  of  Nusselt  number  for  data 
obtained  over  a  wide  range  of  heat  flux,  and  it  was  found  possible  to  corral at a 
thase  values  with  the  kinematic  viscosity  ratio: 


v 

A*  =  30.2  exp  (0.0265  — )  . jg 

b 

After  substituting  this  expression  for  A4  in  the  original  expression  for  eddy 
diffusivity,  it  was  found  that  the  resulting  equation  gave  good  agreement  with 
the  experimental  data  for  l/D  values  greater  than  30,  but  very  poor  agreement 
with  l/D  <  30.  The  latter  case  also  corresponded  to  the  condition  where 
T^  <  Tm  and  it  is  argued  by  Hendricks  and  Hsu  in  their  discussion  of  the 
original  paper  that  this  poor  agreement  can  be  explained  on  the  basis  of  the 
behaviour  of  the  viscosity  ratio,  and  hence  of  the  danping  factor  A4  as 

T^  approaches  T  .  Hess  and  Kunz,  however,  explain  the  lack  of  agreement  as  the 
result  or  an  expended  entrance  region  in  the  heated  tube.  As  with  the  other 
theoretically  derived  relations  for  heat  transfer  to  turbulent  flew  of  super¬ 
critical  fluids,  the  method  involves  complex  and  time  consuming  calculations. 

The  authors  therefore  propose  a  modification  of  the  conventional  equation: 

Nuf  =  F  (fie^,  Pr«)  for  design  purposes,  in  which  the  deviation  between  experi¬ 
mental  and  computed  Nu  values  is  represented  as  a  function  of  the  kinematic 
viscosity  ratio: 


/Equation 
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Nuf  =  0.0208 
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tl  +  0.01457 


The  majority  of  the  data  were  correlated  by  this  relation  to  ±  20}l* 

An  earlier  attempt  to  correlate  data  for  heat  transfer  to  turbulenily 
flowing  supercritical  hydrogen  was  made  by  Hendricks  et  al.  (33)  ana  relies 
upon  a  suggested  similarity  between  forced  convection  film  boiling  and  heat 
transfer  to  fluid  at  near-critical  conditions.  Experimental  data  are  presented 
for  both  these  regions  of  convective  heat  transfer  to  hydrogen  (33*51 )• 

The  scheme  for  correlation  of  supercritical  data  uses  a  pseudo  two-phase 
model  for  the  near-critical  fluid  state.  It  is  postulated  that  the  super¬ 
critical  region  is  an  extension  of  the  two-phase  region,  and  that  fluid  in  the 
near-critical  state  consists  of  an  aggregate  of  two  types  of  molecule  clusters. 
The  fluid  is  taken,  for  simplicity,  to  consist  of:  (a)  a  heavy  species 
corresponding  to  the  tightly  packed  liquid  structure,  with  its  density  (Pme^t) 
equal  to  that  of  liquid  hydrogen  at  the  melting  temperature,  and  (b)  a  light 
species*  with  density  p  corresponding  to  that  of  a  perfect  gas  at  the  bulk 
fluid  temperature,  viz. 


g,b  =  BTh 


The  mass  and  volume  fractions  of  the  respective  species  can  then  be  readily 
obtained  for  supercritical  hydrogen  in  any  condition,  from  a  knowledge  of  the 
density*  a  function  of  the  temperature  and  pressure.  The  bulk  density  is 
given  by 


where  xj  is  the  mass  fraction  of  the  gaseous  species. 

The  mechanism  for  the  boiling-like  process  envisaged  was  described  as 
fellows:  "Instead  of  two  phases,  there  wo'.ild  be  a  continuum  of  densities 
between  a  light  gas-like  species  and  a  heavy  liquid-like  species.  Heavy, 
tightly-packed  clusters  of  molecules  would  migrate  toward  the  wall  and  then 
would  be  broken  up  into  smaller  and  lighter  clusters  and  migrate  toward,  the 
core  of  the  fluid*. 

Hendricks  et  al.  had  been  able  to  correlate  their  data  for  film  boiling 
of  liquid  hydrogen  (51 )  by  means  of  two  parameters:  the  Martinelli  factor, 
hr  “a  a  predicted  Nusselt  number,  Nu’ ,  defined  by 


%  ■  (H-r  C-ST  &rJ 

TPf  Uav  M0'* 

Nu»  =  0.023  [ - - - j  Prf0*4 
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In  Equation  the  fluid  density,  p^,  at  the  film  temperature,  was 
defined'  by 


where  xz  is  the  mass  fraction  of  the  gaseous  species. 

By  plotting  the  ratio  Nu(expt.)/Nu*  as  ordinate  against  the  liartinelli 
parameter  for  the  supercritical  data,  a  very  similar  result  to  vhat  obtained 
with  the  data  for  film  boiling  was  obtained. 

The  correlating  equation  is  of  the  form: 

Nu(expt.)  =  (0.611  7  1.93  Xtt)  . ^ 

The  method  yields  satisfactory  correlation  of  the  experimental  data  over 
wide  ranges  of  heat  flux  and  mass  flow  rate,  although  the  procedure  for 
calculating  heat  transfer  coefficients  is  somewhat  complicated.  The  ratio  of 
experimental  Nusselt  numbers  to  those  predicted  by  the  conventional  "pipe  flow" 
equations  may  have  values  of  up  to  3  in  the  near- critical  region,  giving  support 
to  the  concept  of  an  additional  (boiling-like)  mechanism  enhancing  turbulent 
transfer  of  heat  to  the  fluid.  Use  of  the  pipe  flow  equations  in  design,  while 
these  are  extremely  conservative,  could  lead  to  difficulty  in  the  near- critical 
region  in  that  the  additional  heat  gained  by  the  fluid  from  the  wall  would 
result  in  large  density  changes  and  increase  in  fluid  velocity  in  the  coolant 
passages,  ("vapour  choking")  and  incorrect  estimation  of  the  state  of  the  fuel 
at  the  point  of  injection  prior  to  combustion.  BoWever,  the  conventional  pipe 
flow  equation,  with  fluid  properties  measured  at  "film"  temperature: 

Nuf  =  0.023  Hef°*8  Prf°*'4  . 20 

will  correlate  data  for  near-gaseous  supercritical  hydrogen  to  ±  15£»,  over  wide 
ranges  of  temperature  and  pressure,  provided  that  the  fluid  exhibits  no  unusual 
or  rapid  changes  in  properties  with  temperature. 

An  important  restriction  on  the  pseudo-two-phase  correlation  is  that  the 
bulk  fluid  temperature  oust  be  near  or  above  the  transposed  critical  temperature, 
T  •  In  practice  this  will  not  limit  the  usefulness  of  the  correlation  seriously, 
since  the  critical  temperature  for  hydrogen  is  only  a  little  above  the  normal 
boiling  temperature  (33°  K  and  20°K  respectively)  and  for  those  pressures  where 
the  supercritical  effects  are  likely  to  be  serious,  the  transposed  critical 
temperature  will  be  only  slightly  higher  than  the  critical  temperature. 

A  number  of  correlating  equations  which  have  been  proposed  for  heat  transfer 
to  supercritical  hydrogen  have  been  referred  to  above.  Williamson  and  Bartlit 
(52)  compared  seven  different  correlations  which  hcve  been  proposed  for  hydrogen, 
and  also  the  simple  "pipe  flow"  equation  of  Dittus  and  Boelter,  with  experimental 
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data  for  hydrogen  from  all  known  sources,  covering  a  range  of 'fluid  temperatures 
from  33  to  280°X,  and  pressures  from  1.7  to  97  bars.  They  found  that  although 
individual  correlating  equations  usually  fitted  best  the  data  for  which  they 
were  proposed,  the  best  correlation  for  data  from  all  the  sources  was  the  Dittus 
and  Boelter  equation,  but  with  the  initial  constant  term  halved,  viz. 


Ni^  =  0.012  Re^*8  Prb°*4 


21 


This  correlated  7 Ofi  of  all  the  experimental  data  at  film  temperatures  between 

55  and  280°K,  to  ±  20&  although  scatter  in  the  vicinity  of  the  critical 

temperature  (T  =  33°K)  was  much  greater,  as  would  be  expected, 
c 

For  prediction  of  heat  transfer  coefficients  at  near-critical  conditions, 
where  fluid  bulk  temperatures  are  below  50°K,  it  is  recommended  that  approximate 
values  should  be  estimated  from  available  experimental  data  in  this  region;  or 
from  correlations  such  as  Equation  13,  due  to  Hess  and  Runs,  for  which  some 
reliability  has  been  demonstrated  by  experiment  (38)  at  supercritical  conditions. 

The  Hess  and  Kiunz  correlation  appears  to  give  conservative  predictions  of 
heat  transfer  in  curved  tubes  also,  and  is  recommended  as  a  design  equation  for 
these  geometries  when  modified  to  account  for  the  enhancement  of  mean  heat 
transfer  coefficients  produced  by  tube  curvature. 


5.  SPECIAL  EFFECTS  WITH  HEAT  TRANSFER  TO  LIQUID  HYDROGEN:  PABA-  TO  ORTHO- 
CONVERSION 


The  possibility  exists  that  heat  transfer  to  cryogenic  hydrogen  may  be 
influenced  by  the  conversion  of  the  para-  to  the  ortho-  form  of  hydrogen,  since 
the  equilibrium  composition  of  liquid  hydrogen  varies  with  temperature  as 
follows: 


Temperature,  °K  %  Para-Ha 


%  0rtho-H2 


20  99.8  0.2 

77  50  50 

300  25  75 


Liquid  hydrogen  used  in  operational  rocket  engines  (and  in  the  majority  of  the 
experimental  studies  reported),  is  usually  catalytically  converted  to  give  a 
very  high  proportion  of  para-hydrogen  prior  to  liquefication.  This  prevents 
subsequent  high  evaporation  losses  due  to  the  gradual  conversion  of  the  ortho- 
to  the  para-  form,  with  evolution  of  heat  comparable  to  the  latent  heat  of 
evaporation  of  hydrogen. 

When  liquid  para-hydrogen  at  cryogenic  temperatures  cooes  into  contact 
with  surfaces  at  temperatures  several  hundred  degrees  Kelvin  higher,  rapid 
conversion  of  para-  to  ortho-  could  take  place,  this  representing  an  additional 
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heat  sink.  The  most  probable  mechanism  for  this  conversion  is  a  reaction  at 
the  hot  wall,  possibly  catalysed  by  the  material  of  the  surface,  e.g.  stainless 
steel  or  nickel  alloy. 

Hendricks  et  al.  (51 )  attempted  to  calculate  the  contribution  of  the  para- 
to  ortho-  conversion  to  the  forced  convection  heat  transfer  to  sub-critical 
cryogenic  hydrogen,  assuming  that  the  reaction  was  first  order,  requiring  the 
diffusion  of  the  ortho-  and  para-  species  through  the  laminar  boundary  layer 
at  the  wall.  They  concluded  that  the  conversion  could  have  made  no  perceptable 
difference  to  the  amount  of  heat  transferred  to  the  fluid.  Experimental 
results  with  para-hydrogen  in  which  both  stainless  steel  and  gold-plated  heated 
surfaces  were  employed,  tended  to  confirm  this  finding. 

Thompson  and  Geery  (l5)»  used  lOOjS  para-hydrogen  at  supercritical  pressures 
(?red  and  noted  from  study  of  the  heat  balances  in  the  experiments 

that  the  ratio  ihoreajed  trith  wall  temperature. 

They  concluded  that  an  appreciable  fraction  of  para-hydrogen  was  converted 
to  ortho-hydrogen  at  high  wall  temperatures,  with  absorption  of  the  endothermic 
heat  of  reaction.  The  degree  of  conversion  was  apparently  a  function  of  both 
the  wall  temperature  and  the  residence  time  of  the  fluid  in  the  heated  tube; 
the  latter  was  only  of  the  order  30  milliseconds.  It  is  conceivable  that  heat 
losses  from  the  apparatus  may  account  for  these  observations. 

Perroud  anu  Eebiere  (28)  compared  heat  transfer  to  95/  para-hydrogen  and 
to  60/  para-hydrogen  under  identical  experimental  conditions  and  found  no 
difference  which  could  be  attributed  to  the  para-  to  ortho-  conversion. 

No  3tudie3  are  known  in  which  the  oriho-/para-  composition  of  the  fluid  at 
the  inlet  and  cuclet  from  the  heated  section  was  measured  simultaneously  during 
the  experiments. 


6.  CONCLUSIONS 


It  is  evident  from  this  survey  that  no  consensus  is  obtained  from  the 
published  data  regarding  the  nature  of  heat  transfer  to  hydrogen  under  near 
critical  conditions.  A  number  of  investigations  strongly  suggest  that  under 
these  conditions  the  fluid  behaves  as  if  ir.ro  "pseud-— phases"  exist  during 
turbulent  flow  through  heated  tubes,  with  some  similarities  to  the  behaviour  of 
boiling  liquids.  Additional  evidence  is  available  to  support  this  concept  from 
the  data  for  other  fluids  in  the  same  region  of  reduced  pressure  and  temperature. 
A  convincing  model  for  the  mechanism  of  pseudo-two  phase  heat  transfer  is  not 
30  far  available. 

The  theoretical  and  semi-theoretical  analyses  considered  attempt  to  take 
account  of  the  severe  variations  with  temperature  of  fluid  properties,  and 
discount  the  possibility  of  additional  heat  transfer  from  pseudo-two  phase 
behaviour.  Such  analyses  are  of  limited  practical  usefulness  because  of  their 
complexity  in  application,  and  have  proved  less  successful  in  predicting  heaw 
transfer  to  near-critical  hydrogen  than  to  carbon  dioxide  and  water  in  this 
thermodynamic  region.  The  lack  of  completely  reliable  data  for  a  number  of  the 
thermal  and  transport  properties  of  hydrogen  in  the  critical  region  may  partly 
account  for  this. 
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Untier  most  conditions  of  heat  transfer  to  cryogenic  hydrogen  at  super¬ 
critical  pressures,  the  correlation  of  Hess  and  Kunz  gives  conservative 
estimates  of  convective  heat  transfer  coefficients.  In  general  mean 
coefficients  in  curved  tubes  appear  to  be  up  to  30j£  higher  than  for  the  same 
conditions  in  straight  tubes.  However,  local  heat  transfer  coefficients  at 
the  concave  and  convex  walls  of  curved  tubes  may  be  wid*ly  different,  with 
diminished  heat  transfer  at  the  inside  of  the  curve.  But  it  is  impossible  to 
generalise  these  effects  from  the  available  experimental  data  because  of  the 
influence  of  tubs  diameter  and  other  parameters. 

In  the  case  of  asymmetrical  axial  heat  flux  distribution,  a  severe  minimum 
in  the  local  heat  transfer  may  occur  at  the  position  of  maximum  heat  flux,  and 
this  effect  is  exacerbated  where  the  fluid  temperature  at  this  critical  location 
is  close  to  the  transposed  critical  temperature  of  the  fluid. 

The  available  experimental  data  suggest  that  neither  non-circular  tube 
cro3S-sections,  nor  asymmetry  of  the  peripheral  heat  flux  distribution  produce 
adverse  effects  on  heat  transfer  to  supercritical  hydrogen  when  employed  a3  a 
coolant* 

There  are  no  experimental  data  published  for  heat  transfer  to  hydrogen  in 
multi-tube  assemblies.  With  an  arrangement  of  several  tubes  in  parallel 
between  common  manifolds  flow  instabilities  could  result  for  cryogenic  hydrogen 
at  supercritical  pressures,  in  the  case  of  non-uniform  distribution  between  the 
tubes.  This  has  obvious  application  to  the  cooling  problem  in  a  hydrogen- 
fuelled  rocket  engine,  and  one  would  anticipate  that  the  condition  would  become 
more  severe  as  the  fluid  pressure  approached  the  critical  pressure. 
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NOJfH'ICLATUBI, 


dan;  <  faster  for  grov/th  of  turbulent  eddies 
(Eqn.  It.  a'i  Ref.  58) 

non-dimensional  damping  factor 

specific  heat  at  constant  pressure 

tube  diameter 

effective  diameter  for  non-circular  channel 
geometry 

gravitational  acceleration 

mass  flow  rate 

heat  transfer  coefficient 

enthalpy  of  fluid 

thermal  conductivity 

turbulence  mixing  length  (Eqn.  JO) 

length  along  heated  surface 

constants  in  the  expressions  for 
eddy  diffusivity  (Ref.  17) 

molecular  weight 

pressure 

critical  pressure 

reduced  pressure  (P^^  =  P/P^) 


heat  flux 


generalised  heat  flux 


/  »  •  tf”2\ 


(cal/g  degK) 


(cm/sec2) 

(g/cm2  sec) 
(cal/cm2  sec  degK) 
(cal/g) 

(cal/cm  sec  degK) 


(atm) 

(atm) 


(cal/cm2  sec) 


R^,  R^  upstream  and  downstream  radii  of  curvature  in 
curved  tubes 

R  gas  constant 


(cal/degK  mole) 


temperature 
reduced  temperature 
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critical  temperature 


(°K) 

(°K) 


T 

c 


Tf 

film  temperature  [T^  =  +  Tg)  ] 

(°K) 

Tm 

transposed  critical  temperature  (i.e.  the 
temperature  at  which  the  specific  heat  Cp 
attains  its  maximum  value  at  a  given  supercritical 
pressure) 

(°K) 

<i> 

~x 

reference  temperature  [T^  =  +  x(Tfl  -  T^)] 

for  evaluation  of  fluid  properties  in  variable 
property  heat  transfer  analysis  (fief.  17) 

t- 

(°K) 

u 

velocity  in  stream  parallel  to  heated  surf  ace 

(cin/sec) 

V 

specific  volume 

(cm3/ g) 

XI 

fraction  of  heavy  "liquid*  species  in  near- 
critical  fluid  (Ref.  33) 

X2 

fraction  of  light  "gaseous"  species  (Ref.  33) 

y 

distance  from  heated  surface 

(cm) 

y+ 

P 

K  WJF 1  / 

heat  flux  parameter  =  | '  Cp  -gT~T*~*  J  ( 

'  P« 
Ref. 

V 

It) 

€h 

eddy  diffusivity  (thermal) 

(cm2/sec) 

c 

D 

eddy  diffusivity  (momentum) 

(cm2/sec) 

* 

€ 

m 

momentum  eddy  diffusivity  corrected  :°or  variation 
of  fluid  density  across  stream 

(co2/aec) 

A. 

latent  heat  of  vaporisation 

(cal/ g) 

P 

dynamic  viscosity 

(g/cm  sec) 

P 

density 

(g/co3) 

's 

density  of  perfect  gas  (Ref.  33  and  Eqn.  14) 

(g/ca3) 

*melt 

density  of  liquid  at  the  melting  point 

(g/cm3) 

*>f 

film  density  for  pseudo-two  phase  fluid  (Eqn. 

18) 

(g/cm3) 

r 

shear  stress 

(g/cm  sec2 

V 

kinematic  viscosity  (v  =  p/p) 

(cm2 /sec) 

/Dimensionless 
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Dimensionless  Group 3 


Gr 

No 

Pr 

Re 

Sr 


Grasshof  number  (gr 


X  ~  O  teP, 

'  2  * 
i  I*1 


Nusselt  number 
Prandtl  number 
Reynolds  number 

Sterman  number 


(Nu  =  hD/k) 
(Pr  =  fiCp/k) 
(Re  =  GD/jj) 


Subscripts 
av  average 

b  condition  or  property  in  the  bulk  fluid  stream 

t  film  condition 

S  gas 

s  at  heated  surface 

m  at  transposed  critical  tengierature 


S.  No.  60/69/CJ 


P‘US  **  */m* 


mmn 

mtm 

mmmm 

mumvml 

miiasmiau 

iff«>jafeai 

ii  imMam 

vwts* 


n 

Kljtai 

!»■■■■■ 

laiwmuan 
lawvifeiiH 
iswjfetiH 
muRiiwi 


ii'ttaifctii 


<  UpwiJ 


P -245x10*  N/m* 

1.  O  -  3M  k*/m*  *,  f  -  0-27  MW/m* 
£g-3g  »  •* -o-37  , 

f  0-400  „  ,*-045  * 

*•  »  *t"W*  » 

•  *  » f  —  043  „ 

S.O-281  „  ,f«0* 


Adi  (4)  oad  aoff  (»)  « 
akmgtht  tmbt 


PK3.1.  EFFECT  OF  DIRECTION  OF  FLOW  IN  HEAT 


NE 


VERTICAL  TUBES.FROM  SHITSMAN  (18) 


EXPERIMENTAL  CONDITIONS 


ERPE  1/S/69 


« 

*1: 

«  o 

l|8  6 

JO  ° 

*8s$ 

«  ■  a  • 

0. 


Uj 

5  u. 


fe 


I 

O.  s= 


fc: 

o 

a) 

a* 


* 

o 


g 


*  J  • 

«Q:  b: 


«r  -J 

b  S 

§  §  - 

Kj  £>  *»J 

<  N  <«< 

HI 

i  i 
O I 


hi  o 

S*  a 


3 

ii 

Uj  q: 

u. 

I  I 


s 

I 

Ju 


Q 

2) 

U. 

I 


«*>  *  «o 


FIG. 2.  WAUL  TEMPERATURES  OF  THE  STRAIGHT  PIPES  AT 
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